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Interpolation  with  Curvature  Constraints 


Hafsa  Deddi,  Hazel  Everett,  and  Sylvain  Lazard 


Abstract.  We  address  the  problem  of  controlling  the  curvature  of  a 
Bezier  curve  interpolating  a given  set  of  data.  More  precisely,  given  two 
points  M and  N,  two  directions  u and  v,  and  a constant  k,  we  would 
like  to  find  two  quadratic  Bezier  curves  Ti  and  T 2 joined  with  continuity 
G1,  and  interpolating  the  two  points  M and  N,  such  that  the  tangent 
vectors  at  M and  N have  directions  u and  v respectively,  the  curvature  is 
everywhere  upper  bounded  by  fc,  and  some  evaluating  function,  the  length 
of  the  resulting  curve  for  example,  is  minimized.  In  order  to  solve  this 
problem,  we  first  need  to  determine  the  maximum  curvature  of  quadratic 
Bezier  curves.  This  problem  was  solved  by  Sapidis  and  Prey  in  1992.  Here 
we  present  a simpler  formula  that  has  an  elegant  geometric  interpretation 
in  terms  of  distances  and  areas  determined  by  the  control  points.  We 
then  use  this  formula  to  solve  the  variant  of  the  curvature  control  problem 
in  which  Ti  and  1^  are  joined  with  continuity  C1,  where  the  length  a 
between  the  first  two  control  points  of  Pi  is  equal  to  the  length  between 
the  last  two  control  points  of  ]?2,  and  where  a is  the  evaluating  function 
to  be  minimized. 


§1.  Introduction 

An  important  problem  in  CAGD  is  the  construction  of  curves  interpolating 
given  sets  of  data  that  also  satisfy  constraints  on  their  curvature.  Such  curves 
are  visually  pleasing  and  are  said  to  be  “fair”  [1,2].  Fair  curves  are  also  impor- 
tant in  the  design  of  highways,  railways  and  trajectories  of  mobile  robots  (see 
[9]  and  [6]).  In  these  applications,  curvature  continuous  curves  with  bounded 
curvature  are  desirable.  Constructing  fair  curves  has  been  the  subject  of  recent 
research;  see,  for  example,  [4,5,7]  for  results  about  constraining  the  curvature 
at  the  endpoints,  and  [3,8]  for  results  about  monotonicity  of  curvature. 

In  this  paper  we  consider  the  problem  of  controlling  the  curvature  along 
the  whole  length  of  a Bezier  curve  interpolating  a given  set  of  data.  More 
precisely,  given  two  points  M and  N,  two  directions  u and  v,  and  a constant 
k,  we  want  to  find  two  quadratic  Bezier  curves  T 1 and  T 2 joined  with  continuity 
G1,  and  interpolating  the  two  points  M and  N,  such  that  the  tangent  vectors 
at  M and  N have  directions  u and  v respectively,  the  curvature  is  everywhere 
upper  bounded  by  k,  and  some  evaluating  function,  the  length  of  the  resulting 
curve  for  example,  is  minimized.  We  call  this  the  curvature  control  problem. 
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In  order  to  solve  this  problem,  we  first  need  to  determine  the  maximum 
curvature  of  quadratic  Bezier  curves,  that  is,  to  find  an  exact  formula  in 
terms  of  the  control  points.  Note  that,  for  our  problem,  it  is  not  sufficient  to 
compute  the  maximum  curvature  of  a particular  Bezier  curve  using  numerical 
methods.  Note  also  that  a quadratic  Bezier  curve  is  a parabola  and,  although 
it  presents  no  special  difficulties  to  compute  the  maximum  curvature  of  a 
parabola  in  terms  of  the  coefficients  of  its  implicit  equation,  what  we  require 
is  a formula  in  terms  of  the  control  points. 

In  [8],  Sapidis  and  Frey  give  a formula  for  finding  the  maximum  curvature 
for  quadratic  Bezier  curves.  In  Section  2,  we  recall  these  results  and  present 
a simpler  formula  that  has  an  elegant  geometric  interpretation  in  terms  of 
distances  and  areas  determined  by  the  control  points.  We  then  use  this  formula 
to  solve  variants  of  the  curvature  control  problem.  Definitions  and  motivations 
for  these  variants  are  presented  in  Section  3.1.  We  solve  in  Section  3.2  the 
version  of  the  curvature  control  problem  where  F i and  I 2 are  joined  with 
continuity  C1,  where  the  length  a between  the  two  first  control  points  of  Y 1 
is  equal  to  the  length  between  the  two  last  control  points  of  T2,  and  where  a 
is  the  evaluating  function  to  be  minimized.  In  Section  3.3,  we  prove  that  if  we 
require  in  the  previous  variant  a continuity  G2  instead  of  C1  at  the  junction 
point,  then  there  exist  non-degenerate  data  for  which  there  is  no  solution  to 
the  curvature  control  problem.  However,  if  a solution  exists,  we  show  how  it 
can  be  computed. 

Throughout  the  paper,  curvature  refers  to  non-signed  curvature,  unless 
otherwise  indicated.  We  denote  by  ||pq'||  the  distance  between  points  p and  q, 
and  by  “ x ” and  “ • ” the  outer  and  inner  products,  respectively,  between  two 
vectors. 


§2.  Maximum  Curvature  of  Quadratic  Bezier  Curves 

Let  r be  a quadratic  Bezier  curve  with  control  points  po,  p 1 and  P2  (see 
Figure  1).  Recall  that  T is  defined  for  every  t in  [0, 1]  by  F(t)  = (1  — t)2po  + 
2t(l  — t)pi  + t2p2 . Let  A be  the  area  of  the  control  triangle  P0P1P2  and  m 
be  the  midpoint  of  the  segment  poP2-  We  assume  that  F does  not  degenerate 
into  a line  segment,  i.e.,  po,  pi  and  P2  are  not  collinear. 

Theorem  1.  The  maximum  curvature  of  a quadratic  Bezier  T is  either  equal 
to  ||pim||3/^42  if  pi  lies  strictly  outside  the  two  disks  of  diameter  pom  and 
mp 2,  or  is  equal  to  max-fno,*^}  where  kq  = -^/IIpoPiII3  and  = ^4./ IIP1P2 1| 3 
are  the  curvature  of  T{t)  at  the  endpoints  r(0)  and  r(l). 

Before  proving  Theorem  1,  we  recall  the  result  by  Sapidis  and  Frey  [8] 
characterizing  quadratic  Bezier  curves  with  monotone  curvature. 

Theorem  2 [8].  The  quadratic  Bezier  curve  T has  monotone  curvature  if 
and  only  if  one  of  the  angles  Z(popim)  and  Z(mpiP2)  is  equal  to  or  larger 
than  In  other  words,  F has  monotone  curvature  if  and  only  ifpi  lies  on  or 
inside  one  of  the  two  circles  having  as  diameter  p0m  and  mp 2 (see  Figure  1). 


Interpolation  with  Curvature  Constraints 


193 


pi 


Fig.  1.  The  quadratic  Bezier  curve  T has  non-monotone  curvature  because  pi 
lies  strictly  outside  the  two  circles. 

Sapidis  and  Frey  also  present  in  [8]  the  following  expressions  for  the  max- 
imum curvature  of  quadratic  Bezier  curves.  When  the  curvature  is  not  mono- 
tone along  r,  then  its  maximum  curvature  is  4al/||po-B||3,  where  (see  Figure  1), 
a is  the  distance  between  po  and  p2,  l is  the  distance  between  pi  and  the  line 
joining  po  and  p2,  and  |[p()  B||  is  the  distance  between  p0  and  the  line  passing 
through  P2  and  directed  by  plpo  + plp2-  When  the  curvature  is  monotone 
along  r,  its  maximum  is  reached  at  one  endpoint  po  or  p2  of  the  curve,  and  is 

equal  to  2]]PopiF  or  resPectively- 

We  are  now  ready  to  prove  Theorem  1.  Note  that  the  area  A of  the 
control  triangle  P0P1P2  is  equal  to  alj 2.  Thus,  in  order  to  prove  Theorem  1, 
based  on  the  results  by  Sapidis  and  Frey,  it  suffices  to  prove  that  8-4/||po-B||3  = 
||pim||3/^42  or  2 A = ||pim||.j|poB||.  For  completeness,  we  show  how  our  result 
is  derived  from  Theorem  2. 

We  assume  first  that  pi  lies  strictly  outside  the  two  disks  of  diameter  pom 
and  mp2 . Thus,  the  curvature  n(t),  t € [0, 1],  of  the  quadratic  Bezier  curve  F 
is  not  monotone  by  Theorem  2.  It  follows  that  the  maximum  curvature  of  T 
is  obtained  when  the  derivative  of  n(t)  is  zero. 

The  first  and  second  derivatives  of  the  Bezier  curve  F are 

r'(t)  = 2((1  - f)(pi  - po)  + t(p2  - pi)) 

= 2(pi  - po)  + 2t(p2  - 2pi  + po),  (1) 

r"(t)=2(p2-2p1+p0).  (2) 


The  curvature  of  T at  T(f)  is  thus,  for  any  t € [0, 1], 

= |r'(t)  x r"(t)|  = |4(px  - Po)  x (p2  — Pi)| 

ITOII3  l|r'(t)||3 


giving 


8.4 

rWII3’ 


K{t)  = 


(3) 
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where  A = |(pi  -po)  x (p2  -pi)|/2  is  the  area  of  the  control  triangle  poPiP2- 
The  derivative  of  /«(£)  is 


= — 24.4([|r'(t)||)'  = ~12.4(l|r'(£)||2)' 

K[)  r(£)||4  l|r'(£)||5 

Since  we  assumed  that  the  Bezier  curve  T is  not  degenerate,  p o,  Pi  and  P2  are 
not  collinear  and  thus  A.  ^ 0.  Thus,  k'(£)  — 0 if  and  only  if  ( ||2)/  = 0, 
or  alternatively,  T'(£)  • T"(£)  = 0.  Using  (1)  and  (2),  we  get 

T'(f)  • r"(<)  = 4[(p2  - 2pi  + po)£  + (pi  - Po)]  • [P2  - 2pi  + Po]  = 4(a£  - (3), 

where  a = ||p2  — 2pi  + po||2  and  f3  = — (pi  — Po)  • (P2  — 2pi  + Po)- 

Thus,  the  derivative  of  the  curvature  k(()  vanishes  if  and  only  if  t = r = 
(3 /a.  Note  that  r is  in  (0, 1)  because  the  curvature  of  T is  not  monotone  by 
assumption.  Therefore,  the  maximum  curvature  along  T is  obtained  for  t = t. 

2 A 

Lemma  3.  ||r'(r)||  = 

Proof:  By  (1),  the  square  of  the  first  derivative  of  T(f)  at  r is 

||r'(r)||2  = 4[(p2  - 2pi  + po)r  + (pi  - po)]2  = 4(ar2  - 2 r/3  + ||poPi||2) 

= 4(a^  - 2 £/?  + UpoPiII2)  = -(ollpopill2  - (3 2), 

a a 

where,  as  before,  a = ||p2  — 2pi  +po||2  and  (3  = — (pi  — po)  ■ (P2  — 2pi  + 
po).  Since  P2  — 2pi  + po  = piPo  + P1P2  = 2plfr),  we  get  a - 4||pim||2, 

f3  - — 2poPi  • pirri,  and  thus 

Ilr'(r)||2  = ||Pi^7Z||2 (4IIPirn'll2HP0Pi  II2  - 4(popt  • pirn)2). 

It  follows  from  the  canonical  equation  ( U x V)2  + (U  • V)2  = U~V2,  for  any 
two  vectors  U,  V,  that 


r(r)||2 


4(poPi  x plm)2 

||pim||2 


Now,  |poPi  xpjm|  is  equal  to  A,  the  area  of  the  control  triangle  poPiP2-  Indeed, 
pTrri  = (piPo  +PiP2)/2  and  thus  |p0Pi  x Pim|  = |poPx  x P1P2 1/2  = A.  Thus, 
||r'(r)||2  = 4^42/||pim||2  which  yields  the  result.  □ 

The  expression  of  Kmax  = k (r)  now  follows  easily.  By  Lemma  3,  |]r'(r)||  J 
is  equal  to  8.A3/||pim||3.  Thus,  (3)  gives 

That  ends  the  proof  of  Theorem  1 when  pi  lies  strictly  outside  the  two  disks 
of  diameter  pom  and  mp2 ■ 


Interpolation  with  Curvature  Constraints 


195 


When  pi  lies  inside  one  of  these  disks,  Sapidis  and  Frey  (see  Theorem  2) 
proved  that  the  curvature  of  the  quadratic  Bezier  curve  T is  monotone.  The 
maximum  curvature  is  thus  the  curvature  at  one  endpoint  T(0)  or  r(l).  Equa- 
tion (1)  gives  r'(0)  = 2(pi  — po)  and  r'(l)  = 2(p2— pi).  It  then  follows  from  (3) 


that 


and 


k(1)  = 


A 

lbiP2||3' 


§3.  Controlling  the  Curvature  of  Piecewise  Quadratic  Bezier  Curves 
3.1.  Preliminaries 

Let  r i and  r2  denote  two  quadratic  Bezier  curves  with  control  points  (po,Pi, 
P2)  and  (90,  qi,qi)  respectively,  and  let  F denote  the  concatenation  of  Ti  and 
r2.  The  general  curvature  control  problem  we  address  is: 

Given  two  points  M and  N,  two  unit  vectors  u and  v,  and  a constant 
k,  we  would  like  to  find  two  quadratic  Bezier  curves  F]  and  F2  joined 
with  continuity  G 1 (at  p2  = go ),  interpolating  the  two  points  M and 
N (at  po  and  g2  respectively),  such  that  the  tangent  vectors  at  M and 
N have  directions  u and  v,  respectively,  the  curvature  is  everywhere 
upper  bounded  by  k,  and  some  evaluating  function  is  minimized. 

We  consider  without  loss  of  generality  k — 1;  for  any  k ^ 0,  we  can  obtain  an 
equivalent  problem  where  k = 1 by  scaling  the  plane. 

The  curves  T 1 and  T2  are  connected  (at  p2  = go)  with  continuity  G1  if 
and  only  if  there  exists  p € (0, 1)  such  that  p2  = go  = PPi  + (1  — p)gi.  The 
curve  T interpolates  M and  N,  such  that  the  tangent  vectors  at  M and  N 
have  directions  u and  v,  respectively,  if  and  only  if  po  = M,  g2  = N and  there 
exists  a and  (3  positive  real  numbers  such  that  pi  — po  = ocu  and  g2  — gi  = (3v 
(see  Figure  2).  One  way  to  solve  the  general  curvature  control  problem  is  to 

1)  find  the  set  of  (a,/3,p)  € (0,+oo)2  x (0,1)  on  which  the  curvature  of  T 
is  everywhere  smaller  or  equal  to  1,  and  then, 

2)  find  a value  (a,  /3,  p'j  in  that  set  for  which  the  evaluating  function  is 
minimized. 

In  general,  this  is  a non-linear  optimization  problem  with  non-linear  con- 
straints, and  thus,  cannot  necessarily  be  solved  quickly  and  accurately. 
Clearly,  the  difficulty  depends  on  the  complexity  of  the  set  of  feasible  solutions 
and  on  the  evaluating  function  that  is  to  be  minimized.  Here  we  consider  sim- 
plifying assumptions.  First,  we  require  a continuity  C 1 at  the  junction  point 
between  the  two  curves  T 1 and  T2.  This  fixes  p to  1/2  and  reduces  the  num- 
ber of  variables  to  two.  To  bring  the  number  of  variables  down  to  one,  we 
arbitrarily  consider  a = j3.  We  then  choose  as  evaluating  function  the  length 
a.  By  minimizing  a,  we  ensure  that  all  the  control  points  pi,  p2  = go  and 
gi  remain  close  to  the  the  points  M and  N we  want  to  interpolate;  in  other 
words,  by  minimizing  a,  we  expect  that  the  length  of  the  resulting  curve  F 
will  not  be  too  far  from  its  minimum.  With  these  further  assumptions,  we 
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Fig.  2.  Curvature  control  problem  with  continuity  C1  and  a = ft. 

solve  (in  Section  3.2)  the  given  interpolation  and  minimization  problem,  ex- 
cept for  the  degenerate  case  when  u and  v are  parallel,  for  which  we  prove 
that  a solution  does  not  necessarily  exist. 

In  Section  3.3,  we  also  consider  a = (S,  but  we  require  a continuity  G2 
(instead  of  C1)  at  the  junction  point  between  the  two  curves  T i and  T2.  In 
other  words,  we  require  the  signed  curvature  to  be  continuous  on  T.  The 
variables  are  then  reduced  to  ( a,p ),  but  the  constraint  that  the  continuity  is 
G 2 links  these  two  variables,  and  thus  the  problem  is  actually  one-dimensional. 
We  prove  in  Section  3.3  that  this  set  of  additional  constraints  is  too  restrictive 
in  the  sense  that  there  exists  non-degenerate  data  (M,  N,  u,  v)  that  cannot  be 
interpolated.  However,  if  a solution  exists,  we  show  how  it  can  be  computed. 

3.2.  Curvature  control  problem  with  C1  continuity 

We  consider  here  the  following  variant  of  the  curvature  control  problem: 
Given  two  points  M and  N,  and  two  unit  vectors  u and  v,  we  want 
to  find  two  quadratic  Bezier  curves  T i and  T 2 joined  with  continuity 
C 1 (at  p2  = go ),  interpolating  the  two  points  M and  N (at  po  and  g2 
respectively),  such  that  the  tangent  vectors  at  M and  N have  direc- 
tions u and  v respectively,  the  maximum  curvature  of  the  two  curves 
is  smaller  or  equal  to  l,  the  distances  a = ||poPi||  <md  (3  = ||<7i<?2|| 
are  equal,  and  such  that  a is  minimized. 

See  Figure  2. 

We  show  in  this  section  how  to  solve  this  problem  for  non- degenerate 
data,  that  is  when  u and  v are  not  eollinear.  When  u and  v are  collinear,  we 
show  that  there  is  not  necessarily  a solution. 

As  we  said  in  Section  3.1,  this  problem  is  equivalent  to  finding  the  smallest 
a £ (0,  +oo ) such  that  the  curvature  of  F i and  T2  is  everywhere  smaller  or 
equal  to  1,  where  po  = M,  g2  = N,  pi  = p0  + au,  q\  = q?  — av  and 
P2  = 9o  = (Pi  + ?i)/2. 

We  show  how  we  compute  the  smallest  a € (0,  +oo)  such  that  the  cur- 
vature of  Fx  is  everywhere  smaller  or  equal  to  1.  Computing  the  smallest 
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a € (0,  +oo)  for  1^2  can  be  done  similarly.  We  then  return  the  curve  T defined 
by  the  biggest  of  those  two  a. 

First,  for  any  value  a € (0,  +oo),  we  need  to  determine  an  expression  for 
the  maximum  curvature  of  Tj . By  Theorem  1 , it  remains  to  determine  whether 
the  maximum  curvature  of  Ti  is  given  by  the  maximum  curvature  Kma;c(ri) 
of  the  parabola  supporting  Ti,  or  by  /c0(Fi)  or  the  curvature  of  Fi  at 

its  endpoints  Fi(0)  or  r1(l),  respectively.  Thus,  for  any  value  a € (0, +oo), 
we  want  to  decide  whether  pi  belongs  to  one  of  the  disks  of  diameter  pom 
and  mp2,  where  m is  the  midpoint  of  p0P2  (see  Figure  2).  Let  c and  c'  be  the 
respective  centers  of  these  disks  and  R be  their  radius.  In  order  to  determine 
whether  pi  belongs  to  one  of  these  disks,  we  compute  and  compare  R2  with 
the  distances  ||pic||2  and  ||pic'||2. 

Since  pi  and  qi  are  linear  in  a,  and  p2  = (pi  + gi)/2,  m = (po  + p2)/2, 
c = (po  + m)/ 2,  and  c'  = (m  + P2)/2,  we  have  that  (c  — po)2,  (c  — pi)2 
and  (c'  — pi)2  are  of  degree  2 in  a.  Thus,  R2  < ||pic||2  and  R?  < ||pic'||2 

are  inequalities  of  degree  at  most  2 in  a (namely  a > and 

a2 [(5m  + 3u)2  — (u  — v)2]  — 2a(16u  + 8v)  ■ poq2  + 8||po<72||2  > 0).  By  solving 
these  equations,  we  get  a partition  of  (0,  +oo)  into  two  sets  of  intervals  1 and 
1'  such  that  the  maximum  curvature  of  Ti  is  given  by  Kmai(ri)  for  any  a £ 1, 
and  by  max(/t0(ri),  for  any  a € T . 

With  .A(poPiP2)  denoting  the  area  of  the  control  triangle  P0P1P2,  we  get 
by  Theorem  1,  when  po,  pi  and  P2,  are  not  collinear, 


Kmax 


(ri)2 


IIpiHI6 

•^(P0PlP2)4  ’ 


«0(ri)2  = 


A{p0pip2)2 

IboPill6 


and  Ki(ri)2 


AP0P1P2)2 

IIP1P2II6 


A straightforward  computation  gives 


p\rA  = 


Po<?2  + q(-3m  - v) 


P0P1  = au  and  plpj 


Po<?2  - a(u  + v) 


Thus,  A(poPiP2)  = IpoPi  x P1P2I/2  = \au  x p^q 2 - a2u  x F|/4  and 


/r-  a2  _ («2(3w  + v)2  - 2a(3u  + v)  ■ p^q2  + ||Po92||2)3 
lo(a-iM  x v — au  x poq2) 


Ko(ri)2 


( a2u  xv  — au  x po?!)2 
16a® 


and  «i(ri)2 


4(«2m  Xv  — au  x Poq2)2 
(a(u  + v)-  ptffr)6 


Thus,  Kmol(ri)2  < 1,  «o(ri)2  < 1 and  ki(Fi)2  < 1 reduce  to  inequalities 
in  a of  degree  at  most  8,  6 and  6 respectively.  Finding  the  intervals  of  X 
and  X'  on  which  those  inequalities  are  satisfied  can  therefore  simply  be  done 
by  computing  the  roots  of  the  corresponding  equations.  More  precisely,  the 
smallest  of  (i)  the  smallest  root  of  «max(ri)2  = 1 in  X,  and  (ii)  the  smallest 
root  of  K0(ri)2  = 1 and  «i(Fi)2  = 1 in  X' , is  the  smallest  a for  which  the 
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Fig.  3.  Case  where  po,  p\  and  p2  are  collinear  and  consecutive. 

maximum  curvature  of  Ti  is  smaller  or  equal  to  1.  Such  a solution  exists  when 
5x«  / 0 because  the  maximum  curvature  of  T j goes  from  +oo  to  0 since 
Kmox(r'i)2)  Ko(ri)2  and  «i(ri)2  tend  to  +oo  when  a tends  to  0,  and  tend  to 
0 when  a tends  to  +oo. 

We  have  shown  that,  when  u x v ^ 0,  the  smallest  a 6 (0,  +oo)  such 
that  the  curvature  of  Ti  is  everywhere  smaller  or  equal  to  1,  and  such  that 
the  control  points  po,  pi  and  p2  are  not  collinear,  exists  and  we  can  compute 
it.  Suppose  now  that  there  exists  a £ (0, +oo)  such  that  po,  pi  and  p2  are 
collinear  (see  Figure  3).  Assume  furthermore  that  pi  lies  in  between  po  and 
P2 ; otherwise,  Tj  is  not  smooth  and  does  not  satisfy  the  constraint  on  the 
curvature.  Since  P2  is  the  midpoint  of  pi<?i,  it  follows  that  p0,  pi,  P2  and  q\ 
are,  in  this  order,  on  the  line  L passing  through  po  and  directed  by  u (the 
line  is  necessarily  directed  by  u because  pi  ^ po  belongs  to  that  line).  With 
u X v ^ 0,  qi  does  not  belong  to  L.  Thus,  for  a < a,  the  triangle  P0P1P2  is 
not  flat,  but  tends  to  a flat  triangle,  with  flat  vertex  at  pi,  as  a tends  to  d. 
Therefore,  when  a tends  from  below  to  a,  Tj  tends  to  a straight  line  segment, 
and  the  maximum  curvature  of  Fi  tends  to  0.  Thus,  there  exists  a < a such 
that  the  maximum  curvature  of  Ti  is  smaller  than  1.  It  follows  that  a is  bigger 
than  the  smallest  solution  a we  found  previously.  Therefore,  when  u x v ^ 0, 
there  is  always  an  optimal  solution  with  po,  pi  and  P2  not  all  collinear. 

We  now  show  that,  when  u x v = 0,  there  may  not  exist  a solution. 
Assume  for  example  that  po<?2  is  not  parallel  to  u and  v,  and  that  u + v = 0. 
Then,  when  a tends  to  0,  ko(Fi),  «i(Fi)  and  Kmai(Fi)  tend  respectively  to 
+00,  0 and  +00.  Similarly,  when  a tends  to  +00,  they  tend  respectively  to 
0,  +00  and  +00.  It  follows  that  max(Ko(r1),Ki(ri))  and  Kmax(Fi)  tend  to 
+00  when  a tends  to  0 and  +00.  In  addition,  ko(Fi),  ki(Fi)  and  «;max(Fi) 
are  never  equal  to  0 because  then  ||pim||  = 0 or  A(poPiP2)  = 0 which  would 
imply  that  po,  pi  and  p 2 are  collinear,  which  is  impossible  since  the  two 
rays  starting  at  p 0 and  q-2  with  direction  u and  — v do  not  intersect.  Thus, 
max(Ko(ri),  Kl(ri))  and  Kmax(ri)  are  strictly  greater  than  a positive  constant 
for  any  a € (0,  +00),  and,  by  scaling  the  plane,  this  constant  can  be  scaled  to 
a value  greater  than  1. 


Fig.  4.  Example  where  $ > 0 for  any  a > 0 (po?2  x u > 0,  poqi  x v > 0 and 
it  x v < 0). 
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3.3.  Curvature  control  problem  with  G2  continuity 

We  consider  here  the  following  variant  of  the  curvature  control  problem: 

Given  two  points  M and  N,  and  two  unit  vectors  u and  v,  we  want 
to  find  two  quadratic  Bezier  curves  Tj  and  r2  joined  with  continuity 
G 2 (at  P2  = qo),  interpolating  the  two  points  M and  N (atpo  and  q2 
respectively),  such  that  the  tangent  vectors  at  M and  N have  direc- 
tions u and  v respectively,  the  maximum  curvature  of  the  two  curves 
is  smaller  or  equal  to  1,  the  distances  a — ||poPi||  and  (3  = ||giqf2|| 
are  equal,  and  such  that  a is  minimized. 

As  we  said  in  Section  3.1,  the  problem  is  equivalent  to  finding  the  smallest 
a € (0,  +oo ) such  that  r i and  T2  are  connected  G2  and  their  curvature  is 
everywhere  smaller  or  equal  to  1,  where  po  = M,  q2  = N,  P\  = Po  + au, 
9i  = 92  — av,  and  there  exists  p € (0, 1)  such  that  p2  = q0  = pp1  + (1  — p)qi- 
The  curves  Fi  and  T2  are  connected  G2  if  and  only  if  the  two  signed 
curvatures  of  Ti  and  T2  at  p2  are  equal,  that  is,  by  Theorem  1, 

Pm  X PlP2  _ 9091  X 9192 
2||PiP2||3  2||909i||3  ’ 

when  the  triplets  of  points  (po,Pi,P2)  and  (9o>9i>92)  are  not  collinear.  We 
easily  get  that  plp|  = (1  - p)pl9i,  9o9i  = M l,  Pi9i  - Po92  - a(u  + v), 
PoPi  = OiU  and  q[q)  = av.  Thus,  we  get  that  T is  G2  if  and  only  if 

au  x (1  — p)(po92  - a(u  + v))  _ p(po92  - a(u  + v))  x av 
(1-p)3IIpi9i||3  P3||pi9i||3 

u x po92  — au  x v _ po92  x v — au  x v 
(1-P)2  P2 

p2  - 2/PJ  + $ = 0 where  $ = Po<l2  QU—  (if  p^2  x (u  + v)  £ 0). 

Po92  x(tt  + tO 

Standard  calculations  yield  that  the  equation  p2  — 2 pty  + $ = 0 admits 
a root  in  (0,1)  if  and  only  if  $ € (—1/3,0).  We  can  easily  choose  Po , 92 , w 
and  v such  that  ^ 0 (—1/3,0).  Indeed  (see  Figure  4),  T > 0 for  any  u,  v 
that  are  on  the  same  side  of  Poqi  (i.e.,  po92  x u and  po92  x v have  the  same 
sign)  and  such  that  v lies  in  the  small  wedge  defined  by  poq)  and  u (i.e.,  u x v 
and  po92  x v have  opposite  signs).  We  thus  proved  that  there  is  no  solution 
to  our  curvature  control  problem  for  a set  of  non-degenerate  choices  of  the 
parameters  M,  N,  u and  v. 

However,  when  a solution  exists,  it  can  be  computed  as  in  the  previous 
section.  Indeed,  the  curvature  Kmax(Ti)  can  be  expressed  as  a ratio  of  poly- 
nomials in  a and  p,  and  the  inequality  K.max(Ti)  < 1 reduces  to  a polynomial 
inequality  of  degree  28  in  a.  Similar  remarks  hold  for  kq (T,)  and  /ti(T,). 
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§4.  Concluding  Remarks 

It  remains  open  to  solve  the  curvature  control  problem  when  the  length  of  the 
curve  is  to  be  minimized.  Another  interesting  approach  would  be  to  determine 
how  much  longer  than  optimal  our  curves  are.  Also,  we  would  like  to  consider 
the  case  when  the  data  consist  of  more  than  two  control  points.  Note  also 
that,  because  of  the  high  degree  of  the  equations,  it  is  not  clear  that  the 
solutions  presented  in  Sections  3.2  and  3.3  are  usable  in  an  interactive  curve 
design  context.  This  should  be  tested  with  an  implementation. 
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